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Abstract 

We present a theory for Fano interference in light scattering by individual obstacle, based on a 
temporal coupled-mode formalism. This theory is applicable for obstacles that are much smaller 
than the incident wavelength, or for systems with two-dimensional cylindrical or three-dimensional 
spherical symmetry. We show that for each angle momentum channel, the Fano interference effect 
can be modeled by a simple temporal coupled-mode equation, which provides a line shape formula 
for scattering and absorption cross-section. We validate the analysis with numerical simulations. 
As an application of the theory, we design a structure that exhibits strong absorption and weak 
scattering properties at the same frequency. 
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I. INTRODUCTION 

Understanding light scattering by individual particles is of fundamental importance [lj]. 
In recent year, Fano interference effect 2| in such light scattering has attracted significant 
attention. Individual particles can support resonances. In the vicinity of such resonance, 
the spectrum of either scattering cross-section, or extinction cross-section, can exhibit Fano 
ine shape, where the spectrum varies asymmetrically with respect to resonant frequency 



Such Fano interference effects are also closely related to the classical 



analogue to electromagnetically induced transparency (EIT) [9(, which have recently been 
proposed and demonstrated in scattering systems consisting of metallic elements 1(J. 
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16| . Although the Fano effect in particle scattering has been demonstrated 



in different systems, given the ubiquitous nature of such effect, it would be very useful to 
present a broader theory that captures some of general aspects. 

In this paper, we generalize the temporal coupled-mode theory — previously developed 
to account for Fano interference effect in waveguide or grating systems 17J — to treat light 



scattering by individual resonant obstacle. In the theory of [171], the Fano line shape appears 
in the transmission spectrum of the system. For individual obstacle, transmission spectrum 
is no longer well defined. Instead, the theory here aims to calculate the scattering cross- 
section of the obstacle. Also related to this work, the temporal coupled mode theory has 
recently been generalized to treat small particle scattering in Ref. [3]. The work of Ref. l^ 
however, did not include Fano interference effect. 

The paper is organized as follows: In Section 2, we show that for each angle momentum 
channel, the Fano interference effect can be modeled by a simple temporal coupled-mode 
equation, which provides a line shape formula for scattering and absorption cross-section. 
In Section 3, we compare the theoretical predictions to numerical simulations. Finally, 
in Section 4, as an application of the theory, we design a structure that exhibits strong 
absorption and weak scattering properties at the same frequency. 
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II. THEORY 



A. Brief review of standard scattering theory 

We start by briefly reviewing the standard scattering theory for a single obstacle. For 
simplicity, this paper only discusses the two-dimensional (2D) case where the obstacle is 
uniform in the z direction, but the concept is straightforwardly generalizable to the three- 
dimensional cases. Consider an obstacle located at the origin, surrounding by air. When 
a TM wave (with its magnetic field H polarized along the z-direction) impinges on the 
obstacle, the total field in the air region outside the scatterer can be written as: 

oo 

Htotai = H ° { h t H ?\ k P) exp(iZ0) + KH^(kp) exp(iZ0)) (1) 

l=— oo 

where H is a normalization constant, (p, 9) is the polar coordinates oriented at the origin, k 
is the wave number in air, and H^ (H^) is the /-th order Hankel function of the first (sec- 



ond) kind. Since H\ ' (kp) — > y A e H^-+4) e ^_ — when p — > oo, and taking the convention 
that the field varies in time as exp(— iujt), one can identify h~l and hj as the incoming and 
outgoing wave amplitudes. The power carried in such incoming (or outgoing) wave is 

pt=— \H*m\ 2 ( 2 ) 

(Notice that the power in two-dimension has the unit of W/m.) By choosing the magnetic 
field normalization 

we have then \h~l~\ 2 and \h~[\ 2 representing the incoming and outgoing power carried by 
waves in the /-th angular momentum channel, measured in the unit of W/m. 

In the cases where the scatterer is much smaller than the incident wavelength, or the 
system has cylindrical symmetry, the /-th order incoming wave only excites the same order 
outgoing wave. So we define Ri as 

R, ^ g, (4) 

which can be thought of as a "reflection coefficient" since it relates the outgoing wave to 
the incoming wave in each channel. Moreover, if the obstacle is lossless, the power carried 
by the outgoing wave must be equal to that of the incoming wave. Consequently, 



Ri = e 1 * 1 (5) 



where <pi is a real phase factor. 

We will now calculate the scattering, absorption, and extinction cross-sections of an 
obstacle. For this purpose, imagining a plane wave incident upon the obstacle, the scattering 
(absorption) cross-section is then defined as the total power scattered (absorbed), divided 
by the intensity of the incident field. For 2D system the cross-section has the unit of length. 
Mathematically, for this scenario, we write the total field in the air region outside the 
scatterer as 

Utotai = H I exp(*k -r)+J2 i^H^ikp) exp(z70) J , (6) 

where k is the wave vector of the incident plane wave, Si is related to the scattered field and 
referred to as the scattering coefficient. To connect to Eq. (pQ), the plane wave is expanded 
as 

exp(zk • r) = l \ ( H P( k P) + H ?\ k P)) exp(iZ0). (7) 

l=— oo 

Combining Eqs. <^ and ([7]), and comparing to Eq. (jTJ, we have 

ft = (8) 
Furthermore, in the l-th channel, the scattered power P sct and the absorbed power P a b s are 

Pact. = \S)\ I Hn I 

Following the definition of the scattering cross-section, the contribution to the total scatter- 
ing cross-section from the l-th channel is 

r _ Psct d^N 2 | g o| 2 2A |Q|2 nm 

^SCt,l = -j— - r— — = — Pi| , (1UJ 

In 1 HO \TJ U 71 



and the total scattering cross-section is 



c sct = - V |^| 2 , (ii) 



l=—oo 



where A is the wavelength in air. In the same way, we have the total absorption cross-section 
as 

, oo 

C abs = -— V(Re{^}+|^| 2 ) (12) 

7T z — ' 

/=— oo 
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Finally, the total extinction cross-section, as the sum of the the scattering cross-section and 
the absorption cross-section, is 

2A °° 

C ext = V Re{Si} (13) 

7T z — ' 

(=— 00 

B. Coupled mode theory for light scattering 

The goal of our theory is then to provide a general view about these cross-sections. In 
the following we consider the case that the obstacle supports resonance for the Z-th channel. 
In this circumstance, Fano effect is the result of interference of two pathways: the direct 
"reflectance" of the incoming wave that forms the background, and the outgoing radiation 
from the excited resonance. Using the temporal coupled-mode theory formalism [17 
the dynamic equation for the amplitude c of the resonance is 

% = (-iu - 70 - 7) c + Kh+ 
h~ = Bh + + rjc 

where Uq is the resonant frequency, 70 is the intrinsic loss rate due to for example material 
absorption, 7 is the external leakage rate due to the coupling of the resonance to the outgoing 
wave, B is the background reflection coefficient, and k and 77 correspond to the coupling 
constant between the resonance, and the incoming or outgoing wave, respectively. In this 
section, as we consider only the l-th channel, for notation simplicity we suppress the subscript 
I in all variables in Eq. (114H . Here the amplitude c is normalized such that |c| 2 corresponds 
to the energy inside the resonator [19|]. Note that such coupled- mode formalism is, strictly 



19|, 



speaking, valid only when 70 + 7 <C ujq 



The coupling constants k and rj are related to each other by energy conservation and 
time-reversal symmetry considerations [l^]. In the lossless case (i.e. 70 = 0), if the incoming 
wave is absent (i.e. h + = 0), from Eq. ( Tl4l) we have 

c = Aexr>(— iujnt — jt) 

V ' (15) 

h~ = Ar]exp(—iu t — jt) 

where A is an arbitrary constant. According to energy conservation, the energy leakage rate 
must be equal to the power of the outgoing wave, i.e. 

d|c| 2 



dt 



27|c| 2 = — \h | 2 = — 7777* I c 1 2 (16) 



which requires that 

7/7/* = 2 7 (17) 

Now let us perform a time-reversal transformation for the exponential decay process as 
described by Eq. (Il5p . The time-reversed case is represented by feeding the resonator with 
exponentially growing wave amplitude (h~(—t))*. Such excitation results in a resonant 
amplitude (c(— t))* that also grows exponentially, without the outgoing wave. Described 
this time-reversed scenario using Eq. (fll]) . we have 

Krf = 2 7 (18a) 

Brf + r] = (18b) 
Comparing Eq. fU8al) to (iT7|) . we obtain 

K = Tj (19) 

Eqs. ( TT7|) - (|T9|) are the main results of the Fano interference theory. If the system is 
lossless, the background reflection is B = e 1 ^ and therefore we can determine k and rj as 

K = r] = v ^ e i( * + *- n, ° (20) 

where n is an arbitrary integer. For lossy systems, one should expect that the strongest 
contribution of the loss is to the resonant properties. In such a case, as an ansatz, we 
will introduce a non-zero intrinsic loss rate 70 in Eq. (|14|) . while still approximating the 
background scattering as a lossless process. We will see below in the numerical example 
that such an approach is in fact sufficient for scattering from lossy plasmonic particles. A 
more rigorous approach, which we will not pursue here, will be to treat explicitly the loss 
as coupling to an additional decay channel, the combined system with the additional decay 
channel can then be treated using theory that is similar to Eq. ( 1141) but without the need 
of introducing the intrinsic loss rate. 

We now suppose that the incident wave has a frequency u. The reflection coefficient R 
can be straightforwardly obtained using Eqs. (j!4p and fT20l : 



One can easily identify that in the lossless case (i.e. 70 = 0), the amplitude of R is unity as 
expected. Following Eq. (jSJ), the scattering coefficient is 



s = Ur-i) 



1 (i(^ -^) +70 )( e ^-l)- 7 (l +e ^) 

2 i(oj — w)+7o+7 



(22) 



Substituting Eqs. fl22l) into ffTTj) and ffl2l) . we have the scattering and the absorption cross- 
sections as 



a 



2A 



set 



7T 



1 (i(u 7o )(e^ - 1) - 7(1 + g*) 

2 



o6» 



z(a; - + 70 + 7 

2A 7o7 

7T ( w - u ) 2 + ( 7o + 7) 2 

C. General line shapes in a single channel 



(23a) 
(23b) 



We now use Eq. f l2"3"j) to illustrate general line shapes of scattering and absorption cross- 
section spectra in a single channel. Let us first discuss the scattering cross-section. In the 
lossless case where 70 = 0, Fig. [I] plots the scattering cross-section spectra with different 
phase factor <ft as arised from the background "reflection" . 
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FIG. 1: Scattering cross-section spectrum as given by Eq. (|22j) for the lossless case (70 = 0). (a-c) 
correspond to = 0, tt/2, it, respectively. 



For the phase factor = 0, the spectrum is a Lorentizian function [Fig. 00(a)] . Note 
that the resonance creates maximal scattering at the resonant frequency Uq. Since = 
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FIG. 2: Scattering cross-section spectrum as given by Eq. ([22]) for the lossy case with 70 = 7. 



corresponds to a very weak background scattering, in this case the only contribution to the 
scattering process is from the resonance. An example of a very weak background scattering 
is that for the if -polarization as we consider here, when the radius of a metallic cylinder is far 
smaller than the wavelength. In this case, the presence of a surface-plasmon resonance would 
create a Lorenztian line shape in the scattering cross-section spectrum, and dramatically 
increase the scattering cross-section. 

For all cases where <fi is not equal to or tt, the spectrum is not a Lorentizian, but rather 
exhibits a Fano asymmetric line shape, which exhibits both enhancement and suppression 
of scattering coefficients in the vicinity of the resonance. For example, Fig. [H(b) shows the 
scattering cross-section spectrum for = 7r/2. We note that in the Fano line shape, neither 
the minimum nor the maximum in scattering cross-section occurs at the resonant frequency 
ujq. Through Eq. (123al) . we can determine that the minimum and maximum occur at the 
frequency u>o ± 7tan(0/2) and ojq =f 7cot(0/2) respectively, where the ± sign corresponds 
to < (j) < 7T or —71 < < 0. Therefore, in order to switch the scattering cross-section from 
the minimum to the maximum, the smallest frequency required is 27 which is for the case 
<ft = ±7r/2. This result is in fact consistent with a previous study on optical switching using 



20|. 



Fano resonances 

When = 7r, the background exhibits maximum scattering, and the presence of the 
resonance creates a dip at uq — no scattering occurs at the resonance frequency [Fig. [0(c)] . 
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This effect is closely related to the all-optical analogue of EIT, which occurs when a super- 
radiant state has the same resonant frequency as a sub-radiant state in the same scattering 
channel. The super-radiant state by itself establishes a broad scattering peak and thus can 
be regarded as establishing a background in the vicinity of the sub-radiant resonance. The 
presence of the sub-radiant state introduces a dip and renders the system transparent (i.e. 
with zero scattering) at the resonant frequency. A better treatment of EIT effect however 



2l| 



requires one to treat the super-radiant and the sub-radiant state at an equal footing 
and will be done in later works. 

To show the effect of the intrinsic loss, Fig. [2] shows the scattering cross-section spectra 
when 7q = 7. Comparing it to Fig. dj one can see that the line shapes are similar to the 
corresponding lossless cases, but the presence of the loss reduces the cross-section variation as 
a function of frequency. In particular, the scattering cross section can no longer reach 2\/tt, 
which is the maximal contributions possible in a single channel. Neither can the scattering 
cross-section reaches zero. In the presence of loss it is no longer possible to completely 
eliminate scattering effect at a single frequency through the use of a Fano interference with 
a single resonance. 

We now consider the absorption cross-section spectrum described by Eq. fl23b|) . In con- 
trast to the scattering cross-section spectrum, which varies with the background phase fac- 
tor, the absorption cross-section is independent of the phase factor. Instead, it always has a 
symmetric Lorentzian line shape with its maximum at the resonant frequency. Thus, Fano 
interference effect does not affect the absorption properties of the obstacle. 



III. NUMERICAL VALIDATION 



To check the validity of the above theory, we compare Eq. ( 1231) to the results of numerical 
simulations based on the Lorentz-Mie method of a cylindrical scatterer. Figure [3] shows 
the schematic of the scatterer, which consists of multiple concentric layers of metal and 
dielectric. The permittivity of dielectric is Ed = 12.96, and the metal is described by a 
Drude model e m = l — Up/(u 2 + ijd&) and 7^ are the plasma frequency and the damping 
rate, respectively). In the frequency range of u < u p , a corresponding planar structure can 
support surface-plasmon waveguide modes that are confined in the dielectric region. For 
the cylindrical structure, whispering gallery modes related to such waveguide modes can be 
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FIG. 3: Schematic of a metal-dielectric-metal cylinder obstacle, 
described by a Drude model. The light gray area is dielectric (e^ = 
is assumed as air. 



The dark gray area is a metal 
12.96). The surrounding media 
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formed in the low angular momentum channels 

We start with the lossless case where 7^ = 0. For the geometry parameters of p\ = 
0.285A P , P2 = X p , and P3 = 1.5A P (where \ = 2txc/uj P) with c being the speed of light in 
vacuum), the scattering cross-section for the / = channel is plotted as circles in Fig. H](a). 
It shows a typical Fano resonant line shape around 0.155a;p. 

To compare to the analytic results, we determine the parameters required in Eq. (|23|) 
in the following way: We express the H field in each layer as a linear superposition of 
Jo(kp) and H^\kp), except in the innermost layer where the field is proportional to only 
Jo, and outermost layer where the field is proportional to only Hq . By matching the 
boundary condition at all metal-dielectric interfaces, we obtain a transcendental equation 
with frequency as a variable. By solving for the complex roots of the transcendental equation, 
we determine the resonance frequency (the real part of the root) and the leakage rate (the 
imaginary part of the root). For the structure in shown here, in the lossless case, we obtain 
u>o = 0.1552co> p , 7 = 1.9166 x 10~ 5 u p . The phase factor = — 0.48827T is established by 
calculating the scattering coefficient of a uniform metallic cylinder with the same size. Using 
these parameters, the scattering cross-section calculated by Eq. (I23al) shows an excellent 
agreement with the numerical result [Fig. |4](a)] . 

We now introduce a metal loss characterized by 7^ = 0.001^. Repeating the same 
process as outlined above for the lossless case, we determine a resonant frequency ujq = 
0.1552u;p and a total loss rate 7 + 70 = 1.0492 x 10~ 4 co> p . Assuming that the leakage rate 
is the same as the lossless case, i.e. 7 = 1.9166 x 10~ 5 u p , we have the intrinsic loss rate 
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FIG. 4: Cross-section of the I = channel for the scatterer shown in Fig. O where the geometric 
parameters are p\ = 0.285A P , p2 = A p , and ps = 1.5A P (A p is the wavelength in the vacuum at the 
plasma frequency u p ). (a) Scattering cross-section for the lossless metal case of 7^ = 0; (b-d) the 
scattering, extinction, and absorption cross-section spectra for the lossy metal case of 7^ = 0.001^. 

70 = 8.5749 x lCT 5 ^. The background phase shift <ft is determined using the same procedure 
by considering the scattering of a corresponding uniform metallic cylinder, which results in 
= — 0.48827T + 8.6 x 10~H. Comparing with the lossless case, we notice that the background 
phase shift is not significantly affected by introducing loss in the system. Thus, in the theory 
for the lossy case, we will choose as a parameter the background phase shift <fi = — 0.48 827T, 
and ignore the small loss that occurs in the background. Using these parameters, the 
theoretical spectra again agree well with the numerical spectra [Figs. ID^b-d)]. The example 
here thus provides a validation of the theory in the single channel case. 
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FIG. 5: Scattering (a) and absorption (b) cross-section spectra for the scatterer shown in Fig. [3l 
with the geometric parameters of pi = 0.36A P , p2 = 0.73A P and ps = X p , and the metal damping 
rate of 7<j = 0.001w p . The circles are from the Lorentz-Mie method, and the solid lines are the 
theoretical fit using Eq. (|22j) . 

IV. MULTIPLE-CHANNEL RESONANT OBSTACLE 

The theory developed in Section II can be straightforwardly generalized to include res- 
onances in multiple channels (where the resonance in each channel can only be excited by 
the same channel incoming wave and coupled to the same channel outgoing wave), since the 
total cross-section is the sum of the contributions from all channels. As an example, here 
we apply our theory to a metal-dielectric-metal cylinder scatterer where all channels that 
strongly contribute to the scattering process exhibit the Fano effect. 
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We again consider the structure is shown in Fig. [31 with a new set of parameters p\ = 
0.36A P , P2 = 0.73A P and P3 = X p . The metal has a damping rate 7^ = O.OOltUp. The circles 
in Figs. 0(a), and (b) correspond to the scattering and absorption cross-section spectra, 
respectively, as obtained from numerical simulations. In the frequency range of 0.22 ~ 
0.233a;p, the total cross-section of the system is dominated by the contributions from five 
channels of / = 0, ±1, ±2. For each channel, the scattering cross-section spectrum and the 
absorption cross-section spectrum show the Fano and Lorentizian line shape, respectively. 
We use the above theory to fit these curves. The fitting results are plotted as the solid lines 
in Fig. [5} which again indicates good agreement between theory and simulation. 

We note that this obstacle can be potentially used as a cloaked sensor 23] — when the 
obstacle is placed in an electromagnetic field, it absorbs the energy but creating only mini- 
mum scattering within a narrow range of frequencies. From Fig. [51 one can see that at the 
frequency u = 0.2282co> p , the total scattering cross-section of the obstacle is only 0.03(2A/7r), 
while the total absorption cross-section is 0.32(2A/7r). In the presence of loss, the scattering 
of the particle cannot be completely eliminated in this case, as can be seen from Eq. (12~T1) . 
Nevertheless, the scattering can be substantially reduced as we show here. The Fano in- 
terference effect thus provides an interesting alternative mechanism for creating a cloaked 
sensor as compared to Ref. [231 ] . 



V. SUMMARY AND OUTLOOK 



In summary, we present a theory for Fano interference in light scattering by individ- 
ual obstacle, based on a temporal coupled-mode formalism. We show that for each angle 
momentum channel, the Fano interference effect can be modeled by a simple temporal 
coupled-mode equation, which provides a line shape formula for scattering and absorption 
cross-section. We validate the analysis with numerical simulations. As an application, we 
design a structure that exhibits strong absorption and at the same time weak scattering 
properties. 

We note that this theory is applicable for obstacles that are much smaller than the in- 
cident wavelength, or systems that have two-dimensional cylindrical or three-dimensional 
spherical symmetry. For arbitrary-shaped obstacles or assembled clusters, the overall scat- 
tering matrix of the system is no longer diagonalizable in the angular momentum basis. 
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Instead, for a theory of Fano interference, one may define the channel as the eigenstates of 
the background scattering matrix, the resonance in this case is likely then to couple to several 
scattering channels simultaneously. Such a theory will be developed in future research. 

This work is supported by DARPA/MARCO under the Interconnect Focus Center, by 
AFOSR Grant No. FA9550-04- 1-0437, and by the DOE Grant No. DE-FG 07ER46426. 
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